Conservative parameterization schemes 



Alexander Bihlo^ and George Bluman^ 

^ Centre de recherches mathematiques, Universite de Montreal, CP. 6128, suae. Centre-ville, 
Montreal (QC) H3C 3J7, Canada 

■i- Department of Mathematics, University of British Columbia 
Vancouver (BC) V6T 1Z2, Canada 

E-mail; hihlo@crm.umontreal.ca, ^bluman@math.ubc.ca 

■ Parameterization (closure) schemes in numerical weather and climate prediction models 

■ account for the effects of physical processes that cannot be resolved explicitly by these 
. models. Methods for finding physical parameterization schemes that preserve conservation 

laws of systems of differential equations are introduced. These methods rest on the possibility 
(j^ ', to regard the problem of finding conservative parameterization schemes as a conservation 

^/J ' law classification problem for classes of differential equations. The relevant classification 

problems can be solved using the direct or inverse classification procedures. In the direct 
approach, one starts with a general functional form of the parameterization scheme. Specific 
forms are then found so that corresponding closed equations admit conservation laws. In the 
inverse approach, one seeks parameterization schemes that preserve one or more pre-selected 
O ^ conservation laws of the initial model. The physical interpretation of both classification 

approaches is discussed. Special attention is paid to the problem of finding parameterization 
schemes that preserve both conservation laws and symmetries. All methods are illustrated 

■ by finding conservative and conservative invariant parameterization schemes for systems of 
one-dimensional shallow-water equations. 



^ ■ 1 Introduction 

■ The problem of replacing the continuous governing equations of the atmosphere-ocean system 
by a discrete approximation is that in general no numerical scheme is capable of preserving all 

_ the geometrical features that the initial system of differential equations possesses. Among these 

■ features are symmetries and conservation laws. The lack of a numerical scheme in preserving 

■ fundamental properties of the model has far-reaching consequences on the practical utility of 
the computed results. Simulating the earth system is a relevant but highly complex task and it 
involves an intricate interaction of theoretical insight, data handling and numerical modeling. 

^ ■ Introducing errors in any of these tasks can lead to severe drifts of the forecasted state towards 

^ ! wrong attractors and thus to misleading weather and climate predictions. To remedy this 

challenge for the discretization part of the equations, several structure-preserving numerical 
integrators were developed [17,18,28,29] that might eventually replace the standard integration 
schemes. 

Less research has been carried out so far on the slightly different but related problem of finding 
structure-preserving closure models or parameterization schemes for the subgrid-scale terms that 
inevitably arise when discretizing a nonlinear system of partial differential equations owing to 
the limited resolution one has to employ when integrating a numerical model. Attention to the 
importance of this issue was brought up in [20], where the task of finding invariant subgrid-scale 
closure schemes for the filtered Navier-Stokes equations was investigated. These ideas were 
recently picked up in [6, 24] with the aim of formulating general algorithms for finding local 
parameterization schemes with prescribed invariance characteristics. These methods rely on the 
property that any generic parameterization ansatz, when introduced into an averaged system 
of differential equations turns this system into a class of closed differential equations. There 
exist powerful methods from the field of group analysis of differential equations [4,5,8,10,11,14, 
21,23,25,26] that can be used to study the symmetry properties of such classes of differential 
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equations in an algorithmic way. By a proper interpretation of these methods, one can turn them 
into effective tools for the construction and study of local parameterization schemes preserving 
invariance properties. 

The fact that the parameterization problem in general can be regarded as the study of the 
properties of classes of differential equations is also the key for the analog problem of finding 
parameterization schemes preserving conservation laws. Paralleling the task of finding param- 
eterization schemes that preserve symmetry properties, this problem has important physical 
applications. Preserving physical conservation laws in the parameterization process is a natural 
requirement as there exist various processes that preserve e.g. energy or mass, but that cannot 
be resolved in a particular numerical model and thus have to be modeled in a simplified manner. 
When constructing closure models for such processes it is natural to require the closed system 
of equations to still preserve energy or mass as otherwise the physical consistency of the param- 
eterization scheme is necessarily violated. Developing algorithmic methods that allow one to 
construct such parameterization schemes is thus a worthwhile endeavor. 

A further point in favor of an extension of the tool kit of geometric closure schemes to include 
conservative parameterizations is related to the fact that symmetries and conservation laws of a 
system of differential equations C have different relations to the solutions of £. By definition, a 
symmetry of £ is a property of L itself without regard to posed initial and boundary conditions 
whereas a conservation law is a property holding for every solution of L irregardless of posed 
initial and boundary conditions. In particular, most of the nontrivial symmetries of equations 
in hydrodynamics are broken once classical boundary conditions are imposed (e.g. rigid walls or 
even periodic domains). In contrast, the validity of a conservation law holds in the presence of 
all kinds of initial and boundary conditions. 

Hence it is clear that it is physically relevant to study the relation between conservation 
laws, classes of differential equations and the parameterization problem. Shedding light on this 
relation is the purpose of the present paper. The main result of our study is a step towards the 
first systematic description of general methods for the construction of parameterization schemes 
that preserve or possess particular conservation laws. 

This paper in organized as follows. In Section 2, the theory of conservative parameterization 
schemes is developed. Here, necessary terminology on conservation laws and group classification 
is introduced to demonstrate that the problem of finding conservative parameterization schemes 
can be regarded as a classification problem of conservation laws in classes of differential equations. 
This classification problem can be solved using both direct and inverse methods. These methods 
and the additional requirement on conservative parameterization schemes to be invariant with 
respect to a nontrivial symmetry can successfully be used to narrow down the vast possibility 
one generally has when constructing local subgrid-scale closure schemes for averaged or filtered 
differential equations. As an example, conservative and invariant conservative parameterization 
schemes are constructed for the system of one-dimensional shallow- water equations in Section 3. 
In Section 4, the results of the paper are summarized and suggestions are given for further 
research directions. 

2 Conservative and invariant parameterization schemes 

In this section we introduce some necessary terminology on symmetries and conservation laws, 
which are essential to formulate the theory of invariant and conservative parameterization 
schemes in a proper way. The exposition of the background material follows [1,2,10,21,24-26], to 
which we refer for a more thorough discussion of the underlying notions, methods and theoretical 
concepts. 
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2.1 General notions and statement of the problem 

Let there be given a system of differential equations denoted by £, which consists of L equations 
of the form Ai(a;, u^")) = 0, / = 1,...L, where x = (x^, . . . , rr^*) are the independent variables, 
u = {u^, . . . , u'^) are the dependent variables and n^"^ denote all the derivatives of u with respect 
to X up to order n, with u being included as the zeroth order derivative. 

Definition 1. A local conservation law of the system C \s a divergence expression 

T>j^^\c = (Di$i + • • • + Dp$f)|£ = 0, (1) 

which holds on the solution space of the system C (denoted by l^). The p-tuple of differential 
functions $ = (<I>^(x, n^"*^), . . . , n^™^)), m G No, is a conserved vector of the associated 
conservation law. 

Here and in the following, Dj is the operator of total differentiation with respect to x*, defined 
by Dj = d^i + Uj-du<^, where Uj = d'^'-'^u"' / d{x^y^ ■ ■ ■ d{xPyp, Uj- = duj/dx\ J = (ji, . . . ,jp) is 
a multi-index, ji G Nq and \J\ = j'l + • • • + jp- The summation convention over repeated indices 
is understood. 

Definition 2. A local conservation law of the system C is trivial if the components of the 
conserved vector $ are of the form = {x,u^^^) + i7-'(x, u^"*-*), where the differential 
function vanishes on the solution space of the system C and the differential function is 
a null divergence, i.e. it satisfies ^jH^ = identically. 

As trivial conservation laws satisfy the divergence condition unrestricted by the system of 
differential equations C, they contain no relevant information about the system C. Thus, only 
nontrivial conservation laws are of interest below. 

Definition 3. Two conserved vectors $ and represent the same conservation law (i.e. are 
equivalent) if their difference $ — is a conserved vector associated with a trivial conservation 
law. 

The above definition implies that conservation laws can only be found up to composition 
with trivial conservation laws, i.e. there is not a single canonical representation of one and the 
same conservation law. Thus, formally the space of conservation laws can be defined as the set 
of elements from the factor space of the set of all conserved vectors with respect to the subset 
of trivial conserved vectors. See [26] for more details. 

Conservation laws are conveniently found using the multiplier approach. This method rests 
on an equivalent recasting of the definition of a conservation law (1) in the form 

Dj^^{x,U''"'^) = A'(x, [/ A; (x, (2) 

where the differential functions A = (A"'^(x, f/^^^), . . . , A'^(x, C/*-''^)), r G Nq are conservation law 
(CL) multipliers, also called the characteristics of the conserved vector ^. Note that expres- 
sion (2) holds for arbitrary functions U{x). It is then obvious that for solutions U{x) = u{x) of 
the system C the right-hand side of the above expression vanishes and thus (2) reduces to the 
definition of a conservation law given above, provided that the characteristic A is non-singular 
on the solution manifold of C. 

Expression (2) can be converted into a system of determining equations for the multipliers A. 
This is facilitated by means of the Euler operator. 

Definition 4. The Euler operator with respect to the dependent variable C/* is the differential 
operator given by 

Ei = du. - Bj, djj.^ + D,-, D,,a^.^^,^ - • • • = (-D)-^%. , (3) 
where (-D)-^ = (-Dipi . . . (-Dpp''. 
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The importance of Euler operators in the study of local conservation laws lies in the property 
that they annihilate any divergence expression Dj^-'. In particular the CL multipliers A = 
(Ai(x, t/W), . . . , A^(x, C/W)), r e No, yield a CL oi C if and only of 

Ei(A'A0 = O, i = l,...,q. (4) 

Equation (4) can be split with respect to A; and its differential consequences. This yields an 
over-determined linear system of partial differential equations, which serve as the determining 
equations for the local CL multipliers of the system C. Once these multipliers are found, the 
associated conserved vectors $ can be constructed using e.g. the direct construction method [1, 
2,10]. 

We now move on to the precise statement of the parameterization problem. Let there be 
given a system of differential equations C, Ai{x, u^*^-*) = 0, I = 1, . . . , L and let there be defined 
a filtering operation 

r{u')=u\x)= [ n'{y)G{x,y)dy, (5) 
Jn 

where dy = dy^ ■ ■ ■ dy^ and = / dy. Eq. (5) is the convolution of the variable with the filter 
kernel G = G{x,y). The filter kernel G{x,y) satisfies 



/ G{x,y)dy = l, 
Jn 



see [19,20,27]. This averaging operation can be used to decompose the instantaneous dependent 
variables u according to 

u = u + u' . 

The average u is referred to as the resolved or grid-scale part of the dynamics, while u' includes 
the unresolved subgrid-scale fraction of u. As we do not mix different averaging methodologies 
in one and the same physical problem, we subsequently denote by a bar any mean value of u, 
irrespectively of what averaging operator is used in the concrete problem of interest. 

Example 1. In the classical Reynolds averaging one uses the time average of u, which is de- 
fined by 



1 rh+T 

Vr{u') = u\x*) = lim - / u\t,x*)dt 



where tQ denotes the initial time one starts to average. This time average follows from (5) upon 
setting t = x^, X* = (x^, . . . , x^) and by factorizing 

G{x, y) = G{t-y^)fl G,{x' - y^) = S{x' - x^), 

i=2 i=2 

where Ht is the step function over the interval T and 5 is the delta distribution. The time 
averaging is a Reynolds operator, i.e. it satisfies UiUj = UiUj. Owing to this property, in the 
splitting u = u -\- u' one has u' = since u = u. The average over a product u^u^ thus gives 
u^vJ = u^u^ +u*'m-'', which is the classical Reynolds decomposition that introduces the Reynolds 
stresses u^'u^' into the averaged Navier-Stokes equations. In practical computation a finite 
T < oo has to be chosen and then u = u{t,x*), i.e. the mean value is still time-dependent. 
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Example 2. In large-eddy simulation of turbulence, the classical Reynolds averaging as intro- 
duced in Example 1 is replaced by a spatial filtering approach defined by 

VLEs{u') = u'{t,x*)= [ u'{t,y*)G{x,y)dy, 

for which the filter kernel in (5) is decomposed according to 

p p 
G{x, y) = G{t - y') J] G,{x' - y') = 6{t - y') J] Gi{x' - y'). 

1=2 i=2 

The filters defined in this way are generally not Reynolds operators, i.e. now u' ^ Q as u ^ u and 
thus filtering over products u^u^ produces additional terms of forms not present in the Reynolds 
averaging approach, i.e. u^u^ = u^u^ + u^'u^ + u'^u^' + u^'u^' . 

With the aid of a particular averaging operator (5) the system C is converted into a system for 
the resolved part u, which can be determined by measurements or in the course of a numerical 
simulation of the system C This is done by introducing the splitting u = u + u' into the system 
C followed by an application of a specific filtering (5), which leads to the averaged system of 
differential equations C given by 

AKx,n("),u;) = 0, l = l,...,L. (6) 

In this expression the A;-tuple w = {w^, . . . ,w^) includes all those terms that arise in the course 
of the averaging or filtering and cannot be determined from the knowledge of the mean or filtered 
values u^'^\ For the Reynolds averaging introduced in Example 1, these are terms like u^' vJ' 
(or higher order products as well as their derivatives), while in the case of the spatial filtering 
of Example 2, w would additionally include terms of the form u^'& , etc. This is essentially 
the closure problem, i.e. Eqs. (6) include more unknown than known quantities. Thus, as they 
stand Eqs. (6) cannot be used for a numerical integration unless one expresses the additional 
unknowns w in terms of certain known expressions. 

Definition 5. A local parameterization or subgrid-scale closure model assumes a functional 
relation between the unknown subgrid-scale terms w and the mean values u^^'\ r S Nq i.e. 

= f{x,u^'-^), i = l,...,k, (7) 

for certain parameterization functions f = (/^, . . . , f''). 

Introducing a local parameterization scheme (7) into system C (6) leads to a closed system of 
differential equations for the mean values u. The inherent problem of this construction is that 
in most cases of interest the information contained in u and its derivatives is not sufficient to 
determine the entire subgrid-scale structure contained in w. The art of constructing physical 
parameterization schemes is to determine the parameterization functions / in such a manner 
that the assumption (7) will allow one to find u from the closed system C with sufficient accuracy. 

Finding suitable parameterization functions / that lead to realistic results for u can be 
rather tedious. One general methodology to restrict the vast number of possible forms for the 
parameterization functions is to choose them in such a manner that the resulting closed system 
of differential equations preserves certain nontrivial geometric properties such as conservation 
laws and/or symmetries. This motivates the following definition. 

Definition 6. A local parameterization scheme is called conservative provided that the closed 
class of differential equations C preserves certain nontrivial local conservation laws. A local pa- 
rameterization scheme is called invariant provided that the closed class of differential equations 
C preserves a nontrivial point symmetry group G. 
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Conservative parameterization schemes can be found using techniques analogous to those for 
the group classification of classes of differential equations. A class of differential equations C\s 
is a system of differential equations of the form Ai{x,u^'^\0(x,u^'^^)) = 0, I = 1, . . . , L, which 
is parameterized by a /c-tuple 9 = {6^, . . . , 6^) of differential functions that satisfy a system of 
K £ No auxiliary differential equations of the form Sj{x,u^'^\9^"^\x,u^"'^)) = 0, j = 1, . . . , K , 
the solution set of which is denoted by S. The system of auxiliary equations in part specifies 
the properties of the class and it is regarded as a system for 9, i.e. x and u^"^^ play the role of 
independent variables. To complete the description of the class C\s one usually takes into account 
a constitutive inequality, S(x, n^"), 0(™)(x, n^"))) ^ 0, which guarantees that all equations from 
the class share some joint properties (e.g. a particular derivative does not vanish, all equations 
of C\s are linear or nonlinear, etc.). 

It is the purpose of conservation law classification to systematically investigate the CLs of a 
class of differential equations. By substituting the general closure scheme (7) into the averaged 
system £ (6) one obtains a class of closed differential equations for u, 

A,(x,n("),/(x,nW)) = 0, l = l,...,L, (8) 

in which the parameterization functions / play the role of the arbitrary elements 9. In the 
following subsections we will introduce and discuss methods that allow one to specify the pa- 
rameterization functions / in such a manner that the system (8) has certain nontrivial local CLs. 
The corresponding classification methods for finding parameterization schemes that possess non- 
trivial maximal Lie invariance groups where introduced in [11], see also [24]. Combinations of 
invariant and conservative parameterization schemes are also possible and will be discussed at 
the end of this section. 

Remark 1. A problem related to the search for parameterization or closure schemes for the 
subgrid-scale terms arising in averaged differential equations is to search for extensions of differ- 
ential equations that preserve some of the geometric features of the original differential equations 
C. Physically, such extensions could be e.g. adding dissipation terms to a non-dissipative model 
or source or sink terms to transport equations. Depending on the nature of the included pro- 
cess, the addition of such extra terms may alter the structure of the initial system of differential 
equation but might still retain some of the geometric features of the original model. A possible 
research question is thus to construct extra terms for the system C: A;(x, n^"^) = 0, I = 1, . . . L 
in such a manner as to preserve certain CLs and/or symmetries of C Mathematically, this 
is done by investigating systems of the form A;(a;,n^"^) = gi{x,u^^^), for a certain L-tuple 
g = {g^{x, u^^^), . . . g^{x, u^^^)) of extensions. It is obvious that this system can be brought into 
the form (8) if the functions / are interpreted as the additional terms that extend the initial 
system £ and no averaging operation is involved, i.e. u^^^ = u^") . Consequently, the same meth- 
ods as introduced below for solving the parameterization problem for system (8) can be used to 
tackle this kind of problem. 

2.2 Conservative parameterizations via direct classification 

In order to discuss the method for finding conservative parameterization schemes, the following 
definitions are useful. 

Definition 7. An equivalence transformation ip from the class C\s is a point transformation 
on the space {x,u^'^\9), which is projectable on the spaces of (x, < n' < n, such that 
\/9 £ S : 9' = f9 £ S and the restriction of ip to the space of (x,u^"^), denoted by '/'l(x,?i("))' a 
point transformation from Cg to Cql Here, Cg and Cg' are equations from the class C\s. 

Thus, equivalence transformations are point transformations that map one system of dif- 
ferential equations from a given class C\s to another system of differential equations from the 
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same class. The collection of all equivalence transformations forms a group, which is called the 
equivalence group C. 

Definition 8. Let there be given two systems of differential equations from the class C\s, 
denoted by and Cg', which have CLs with conserved vectors $ and respectively. The 
pairs (Cg,^) and {Cg',^') are equivalent with respect to the equivalence group C if there 
exists a point transformation if £ G~ that transforms the system Cg to the system Cg/ and 
which transforms the conserved vector $ in such a manner that ^ = (p{x,u^^\^) and are 
equivalent conservation laws, see Definition 3. 

In this definition, the action of a point transformation ip 6 G~ on a conserved vector ^> has 
the explicit form 

¥{x,u^''^) = (a;, i = l,...,p, 

\DxX\ 

where |Da;2;| is the determinant of the matrix (D^jx*). See [9,10,26] for more details. 

The direct classification procedure for finding parameterization schemes with prescribed CLs 
can be formulated in the following way. For a given fixed general form of the parameterization 
functions /, determine those CLs that hold for any equation from the class (8) (i.e. for all 
admissible forms of /) and find all the inequivalent equations from that class that have additional 
CLs. 

To make the classification problem tractable, one first chooses the general form of param- 
eterization functions / one aims to study, i.e. one determines which variables x and u^"^^ the 
functions / should depend on. This choice is physically motivated. Once the general form of 
/ (hence the system of auxiliary equations S) is fixed, one can solve the classification problem 
taking into account the equivalence of CLs as embodied in Definition 8. This means that one 
determines the equivalence group of the general class of closed differential equations of interest 
and then solves the determining equations (4) for CL multipliers. One seeks those values of 
/ (up to equivalence) for which the determining equations for CL multipliers yield additional 
multipliers beyond those for generic parameterization functions. An example for this procedure 
is given in Section 3.1. 

The direct group classification method yields a list of inequivalent equations C j from the class 
C\s that possess inequivalent nontrivial local CLs. Using this list of all possible conservative 
parameterization schemes from the predefined class, one can then test the different schemes 
obtained and select the most appropriate one as a candidate closure scheme for the process of 
interest that needs to be parameterized. 

Physically, the method of finding conservative parameterization schemes using the direct 
classification approach might be most appropriate in the case when one seeks to represent 
processes that are not already included in the dynamics resulting from the system C The 
reason for this is that by means of the direct classification method one might obtain closed 
differential equations that have CLs not possessed by the original system C. 

2.3 Conservative parameterizations via inverse classification 

A different ideology for finding parameterization schemes is the following. Assume that the 
original system of differential equations C has a certain number of nontrivial local CLs. The 
averaging of a differential equation certainly disturbs the geometric structure of the equation 
but it might nevertheless be desirable that the averaged system share some CLs of the original 
system of equations. An example for this is a process that conserves energy but needs to be 
parameterized in a given system of differential equations. For the sake of physical consistency, 
the closed differential equations should conserve energy and thus only parameterization schemes 
that are compatible with energy conservation can be considered. 
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This discussion is related to what is called the inverse classification problem and in the 
framework of a conservative parameterization scheme, it can be formulated in the following way. 
Let there be given an initial system of differential equations C. One first determines CLs holding 
for the original system of differential equations C: A/(a;,n(")) = 0, / = 1, . . . L through the CL 
multiplier approach. Depending on the complexity of the problem of interest, one might not be 
able to obtain an exhaustive description of all CLs but rather restricts oneself to CLs associated 
with characteristics A{x,u^^^) for some fixed (often low-dimensional) r. Among the CLs of C, 
one selects, using physical reasoning, the associated multipliers of those CLs that one aims to 
preserve in the course of the parameterization process. As in the case of the direct classification 
method, one next averages the system C and determines the general functional form of the 
parameterization functions / to be used in the class of parameterization schemes (7). The 
final step is to plug the class of averaged closed differential equations (8) into the determining 
equations (4) for CL multipliers. Since the multipliers that the resulting equations from the 
class (8) should admit are fixed, the determining equations for local CL multipliers are thus 
converted into a system of determining equations for the parameterization functions /. Solving 
this system leads to all equations from the class £.\s that have the same CL multipliers A(x, u^^'^) 
as the original system C, with u^^^ being replaced by the corresponding averaged values u^^\ 

A nontrivial question in this construction procedure is to determine in advance whether at 
least some systems from the class C\s selected has the CLs associated with the chosen multipliers 
A(x, u^^^) stemming from the original system C, i.e. whether the determining equations (4) yield 
any nontrivial solutions. A natural strategy to overcome this problem of possibly triviality of 
the solution of (4) is to (i) either assume that the class of closed equations C\s is rather wide 
(i.e. that the function / depends on a large number of variables from u^^'^) or (ii) to only require 
a suitable small set of CLs being preserved by the parameterization scheme. From the physical 
point of view, the first strategy should be the method of choice. 

Although this possible triviality of the solution of the determining equations for the parame- 
terization functions and the associated failure in finding nontrivial conservative parameterization 
schemes seems undesirable, it nevertheless includes important physical evidence. It simply means 
that for the parameterization ansatz selected, no element of the class of closed equations can 
satisfy the requirement of retaining the desired CLs and thus might indicate that the initial 
parameterization ansatz was flawed. 

In contrast to the direct classification method, which might lead to conservative parame- 
terization schemes that possess CLs which do not hold for the original system C, conservative 
parameterization schemes derived using the inverse classification procedure by construction yield 
no additional CLs for the resulting closed system of differential equations. The inverse classi- 
fication strategy might thus be best suited for processes already included in the full dynamics 
of the original system £, but that cannot be explicitly resolved because of e.g. computational 
limitations. 

Remark 2. The existence of at least the trivial solution of Eqs. (4) in the case when the CL 
multipliers are fixed and possible forms of / are sought, i.e. that the system of determining 
equations is compatible is guaranteed by the fact that (8) can be rewritten in the form 

A,(a;,n('^)) = 5z(x,n«), l = l,...,L, 

where the left-hand side is the same as in the original system £ provided that is used in place 
of u^'^^ and the right-hand side is a functional combination of the parameterization functions /. 
Thus if 5 = (and hence / = 0), the above equation admits the same CL multipliers A as for 
the system C in which is used instead of u^^' . 
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2.4 Conservative and invariant parameterizations 

Methods for finding parameterization schemes with symmetry properties using methods from the 
group analysis of differential equations were introduced in [6,11,24]. There is neither a practical 
nor a theoretical objection against a parameterization scheme preserving both invariance and 
CL properties. Indeed, the compatibility of these two concepts was explicitly demonstrated 
by constructing CL and invariance preserving hyperdiffusion schemes for the two-dimensional 
barotropic vorticity equation on the beta-plane [6]. We now outline how to systematically 
construct invariance and CL preserving parameterization schemes. 

As in the case of conservative parameterization schemes, two main methods are applicable 
to determine parameterization schemes with symmetry properties. These methods are straight- 
forward applications of the group analysis of differential equations and the required techniques 
are based on the direct and the inverse approach to the symmetry classification problem, re- 
spectively. The key to the construction of conservative invariant parameterization schemes is 
that the closure models resulting from the conservative parameterization procedure as outlined 
in the previous two subsections are usually still classes of differential equations. These classes 
are generally narrower than the initial class given by Eqs. (8) but nevertheless include arbitrary 
constants or parameter functions with respect to which the usual symmetry classification prob- 
lem can be carried out. We only outline the main ideas of this construction below as a more 
detailed exposition of the methods available in the field would require a substantial extension of 
the text. Moreover, the symmetry analysis and group classification of differential equations is a 
well-investigated subject. See [1-5,8,10,11,14,21-23,25,26] and references therein. 

Special attention will be also paid to models that are derivable from variational principles, 
i.e. which are Euler-Lagrange equations. For such systems there is a close connection between 
symmetries and CLs that can be utilized to construct invariant conservative parameterization 
schemes. In particular, in this situation, all CL multipliers are symmetries but the converse is 
false [1,2,8,10,21]. 

Invariant conservative parameterizations via direct symmetry classification. In the 

direct symmetry classification approach one starts with a given class of differential equations 
and determines those symmetries that hold for all equations in the class. These symmetries form 
the kernel of maximal Lie invariance groups of equations in the class. The direct classification 
problem is solved by finding all equations in the class that have symmetry extensions with 
respect to the kernel and this investigation is carried out up to the equivalence imposed by the 
equivalence group 

Depending on the complexity of the class of closed differential equations (8) with conservative 
properties, different strategies for solving the direct symmetry classification problem can be 
adopted. For simple classes depending only on a few arbitrary constants or parameter functions 
with few arguments, the direct integration of the determining equations of Lie symmetries up to 
equivalence using compatibility analysis is the method of choice. This method yields a complete 
list of inequivalent equations from classes of the form (8) that are both conservative and have 
nontrivial symmetry properties. If the structure of the class of (8) is too complicated for a 
direct integration of the determining equations of Lie symmetries then the algebraic method of 
group classification can be used. With this method one aims to find symmetry extensions of 
the kernel that are induced by transformations from the equivalence group of the class under 
consideration. The algebraic method thus reduces the problem of finding symmetry extensions 
to the problem of finding inequivalent subgroups of the equivalence group Similar to the 
direct integration of the determining equations, the algebraic method of group classification can 
lead to the complete solution of the group classification problem, namely for classes of differential 
equations possessing the normalization property, see [25]. If the given class of equations is not 
normalized, then the algebraic method will not lead to a complete description of all possible 
inequivalent symmetry extensions of the kernel. One will still find those symmetry extensions 
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that are induced by the equivalence transformations of the class of equations of interest but there 
can be other symmetry extensions of the kernel that cannot be found from the classification of 
the subgroups of the equivalence group. The algebraic method of direct group classification 
for classes that are not normalized is also known as preliminary group classification. A more 
detailed discussion of the techniques available in the field of direct group classification can be 
found in [5,14,25]. 

Irrespective of what method is used to (partially) solve the direct symmetry classification 
problem, all the systems of closed differential equations obtained in the classification procedure 
have the same CLs but different (inequivalent) maximal Lie invariance groups. The resulting 
closure schemes then have to be tested numerically to find the most suitable invariant and 
conservative representation for a given subgrid-scale process. 

Invariant conservative parameterizations via inverse symmetry classification. The 

inverse symmetry classification problem is to find all those equations that have a prescribed 
symmetry property. See [11] for the first systematic outline of this problem for both ordinary 
and partial differential equations. Here, rather than starting from a given class of differential 
equations and describing the invariance properties of equations from this class as done in the 
direct symmetry classification, in the inverse classification one starts with a given Lie group of 
transformations and seeks to find the class of equations (up to some order n) that is invariant 
under the selected group. The inverse symmetry classification procedure rests on the following 
theorem [8,10,11,21,22]: 

Theorem 1. Let there be given a Lie group of transformations G acting on a manifold M . 
If the nth prolongation of G acts regularly on the nth order jet space J" and if there exists a 
functionally independent system of nth order differential invariants h, ■ ■ ■ ,In, h = Ii{x,u^^^), 
then any nth order system of differential equations C admitting G as a symmetry group can 
be rewritten in terms of these differential invariants, i.e. Ai{x,u^'^^) = Ai{Ii, . . . , = 0, 



This result is known as the replacement theorem [12,16]. It implies the existence of certain 



One can use the replacement theorem to construct parameterization schemes that have spec- 
ified CL and invariance properties. To this end, one first determines the complete system of sth 
order differential invariants li of the maximal Lie invariance group G of the original unaveraged 
system of differential equations C or of an appropriate subgroup G^ C G, where s = min(n, r), 
n and r being the highest order derivatives arising in £ and the parameterization ansatz (7), 
respectively. These invariants can be found either using infinitesimal techniques [8,10,11,21,23] 
or moving frames [12, 15, 16]. 

Secondly, once the invariants of G (or G^) are known one can find the multipliers Tf and 
thus obtain the invariant representation (9) of the system C Now suppose that the parameter- 
ization functions / in the class (8) have been determined in such a manner that the resulting 
equations from (8) have the desired CLs. As in Remark 2, one represents (8) in the solved form 
Ai{x, u^")) = gi{x, u^^'^), I = 1, . . . , L. These closed equations will be invariant under G (or G^) 
for those functions g that satisfy the system of equations 



1 = 1,. ..,L. 




rf{x, u("))g.(x, n('-)) = ~gi{h, ...,In), 1 = 1,... 



L 



(10) 



using the same multipliers Tf (replacing u' 
differential invariants /j. 



with M*-"^), for some A^-tuple of functions g of the 
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Similar to the discussion in Section 2.3, it is a nontrivial question to determine in advance 
whether condition (10) can be satisfied for a given set of conservative parameterization functions 
/ and a chosen symmetry group G (or G^). In general, the wider the symmetry group G (or 
G^) is, the more specific are the forms of the differential invariants li and hence the more 
general the class of conservative parameterization schemes has to be in order to allow jointly 
for CL and invariance properties. On the other hand, if from the physical point of view it is 
required that the parameterization of a certain process should have specified invariance and 
CL properties, then failure in satisfying condition (10) may again point to an inappropriately 
chosen parameterization ansatz for / and thus is one more check for the consistency of a physical 
parameterization scheme. 

Variational parameterizations for Lagrangian systems. As mentioned previously, there 
is a direct close connection between symmetries and CLs for systems of differential equations 
that can be derived from a variational principle. More precisely, for each one-parameter Lie 
group of point transformations (or, more generally, one-parameter group of higher order local 
transformations) that leaves invariant the functional 

S[u] = [ L(x,n("))dx, (11) 
Jn 

where dx = dxi ■ ■ ■ dxp, to within a divergence, there is a local CL of the Euler-Lagrange equa- 
tions associated with (11). This result is the celebrated Noether theorem [8,10,21]. Practically, 
the invariance of the functional S[u\ under point transformations of the form x = x{x,u,e) and 
u = u{x,u,e), with associated infinitesimal generator Q = ^*(x,n)9^.i + (f)°'{x,u)dua , can be 
determined by checking whether the condition 

Q(")(L) + LD,f = DiK' (12) 

is satisfied, where Q*-"^ denotes the nth prolongation of Q and B = {B^(x, n^™)), . . . , B^{x, u^™"^)) 
is some p-tuple of differential functions. In this case, it can be proved [8, 10, 21] that the 
characteristic rj = (t/^, . . . ,7]'^) of the vector field Q written in evolutionary form, which is given 
by T]" = (j)"^ — £,^uf is also a set of q local CL multipliers of the associated Euler-Lagrange 
equations of 5[n], i.e. ij = A. This means that a set of local conservation CL multipliers satisfies 
the determining equations for local symmetries, in evolutionary form, of the corresponding 
Euler-Lagrange equations. 

In practice, the Euler-Lagrange equations follow from taking the variational derivative of (11), 
which boils down to applying the Euler operator E„ to the Lagrangian function L. The follow- 
ing definitions and theorem will be useful for the construction of variational parameterization 
schemes, i.e. parameterization schemes that preserve the variational structure of a system of par- 
tial differential equations. A more comprehensive discussion on the relation between symmetries, 
CLs and variational forms can be found in [8,10,21]. 

Definition 9. Let there be given a system of differential equations C: Ai(x,u^"'^) = 0, I = 
1,...,L. The linearizing operator {Frechet derivative) associated with C is the matrix-valued 
differential operator Dc whose components are given by 

Definition 10. The adjoint of the linearizing operator D£ is the matrix-valued differential 
operator whose components are given by 

{DIUV^ = (-D), [^V^) ■ 

for any differential function V = (V'^ , . . . ,V^), where = V\x,u^'^'^). The operator D£ is 
self-adjoint if and only if = D^. 
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Theorem 2. Let C: A;(x,u("')) = 0, Z = 1, . . . , L 6e a system of differential equations. The 
system C can he derived from a variational principle (11) if and only if the linearization operator 
Dc associated with C is self-adjoint, i.e. D*^ = Dc- 

Theorem 2 is the key for the construction of variational parameterization schemes. Suppose 
that the system C admits a variational form (11). Averaging of the equations from system C 
and assuming a general parameterization ansatz (7) leads to a system of equations that can be 
brought into the form 

A,(x,nW) = <7,(x,nM), / = 1,...,L. (13) 

See also Remark 2. Since the linearization operator associated with the left-hand side of the 
above expression is by supposition self-adjoint, the above closed system of differential equa- 
tions (13) will remain variational if and only if 

D* = D„ (14) 

i.e. the linearization operator associated with the right-hand sides of system (13) must be self- 
adjoint. This imposes the required conditions on the function g = [g^ , . . . ,g^) to yield a varia- 
tional parameterization scheme. 

Since the one-dimensional shallow- water equations (17), which will be our running example 
for finding conservative parameterization schemes, are not derivable from a variational principle, 
we illustrate the idea of variational parameterization schemes for the potential Korteweg-de Vries 
equation. 

Example 3. The potential Korteweg-de Vries (KdV) equation is obtained from the usual KdV 
equation vt + vvx + Vxxx = by the differential substitution v = Ux, i.e. 

\Kdv{x, M^^)) = Utx + UxUxx + Uxxxx = 0. (15) 

This system is variational with L = ^u^x ~ g^x ~ \utUx being the Lagrangian. The maximal 
set of point symmetries of (15) is infinite dimensional and spanned by the vector fields 

dt, dx, tdx + xdu, j{t)du, 3tdt + xdx - udu, 

where 7 = 7(t) runs through the set of smooth functions of t. The first four vector fields satisfy 
the variational symmetry condition (12) and are thus associated with local parameterization 
laws. The characteristics of these symmetries are (i) rj = —ut, (ii) i] = —Ux, (iii) t] = x — tUx 
and (iv) rj = j{t) to which the following CLs A(x, U, Ux, Uxx, Uxxx)^pKdv{x, U^^^) = Dfp + DxF 
correspond [10]: 

A p F 



Ux ~\~ '^Uxxx 


ir/S 1 TT TT 


4 "-^x '-^x '-'xxx 


UxUtxx ~\~ Uxxx ~^ UxxUtx, 


Ux 




In^ 4-TT TT — 
3^x ' '^x'^xxx 


2 xxi 


X — tUx 


\tUl + xUx 


■^tUx tUxUxxx ~l~ 2t^xx ~l~ 2"^^^ Uxx ~\~ xUxxx: 


lit) 


lUx 


\-iUl + -iUxxx - 


ItU. 



We Reynolds average the potential KdV equation to obtain utx + UxUxx + Ux 



To keep the computations as simple as possible we aim at finding the unclosed term 
function of Ux and Uxx in such a manner that the closed equations from the class 

Utx ~\~ UxUxx ~\~ Uxxxx — diu, Ux, Uxx), 
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are still derivable from a variational principle. Prom now on we omit bars over u to simplify 
the notation. The linearization operator associated with the right-hand side of this closed 
class of potential KdV equations is given by D = g(„ + Qu^D^ + S'li^a.-D^i which has the adjoint 
D* = 5n,,Dx + i^^xgu^^ -gujB^ + Dlgu^^-B^gu^ +gu. Thus, condition (14) yields the system 

Dx5«.. - gu^ = 0, D^5„^^ - B^gu, = 

that must be satisfied by function g for an equation from the closed class of potential KdV 
equations to remain variational. Solving this system of differential equations yields g = {c2Ux + 
CijUxx + g^{u)^ where Ci and C2 are arbitrary constants and g^ is an arbitrary function of u. 
Hence 

Utx + UxUxx + Uxxxx - {c2Ux + Ci)Uxx " 5^ = 0, (16) 

is the only admissible form of equations from the class of closed potential KdV equations that 
is an Euler-Lagrange equation. The associated Lagrangian of (16) is given by 

where G{u) = J g^du. For arbitrary g^, the only point symmetries admitted by Eq. (16) are 
generated by the vector fields dt and dx, which are variational and give rise to two local CLs 
of (16). 



3 Conservative closure schemes for the shallow-water equations 

In this section we construct conservative parameterization schemes for the one-dimensional sys- 
tem of shallow-water equations. In nondimensional form, this system reads as 

Ai{t,x,u^'^\h^^'>) = ut + uux + hx = 0, ^^^^ 
A2{t,x,u^^\h'^^^) = ht + uhx + hux = 0, 

where u is the velocity and h is the height of the water column over a fixed reference level. For 
the sake of simplicity, we assume a flat bottom topography in which case h can be considered as 
the total height of the water column. The first equation is the momentum equation, the second 
equation the shallow-water continuity equation. 

Reynolds-averaging the above system (17) with the averaging interval T being finite and 
following the averaging rule for products, u^u^ = u^u^ + u^'u^' , then leads to the averaged 
shallow-water equations 

ut + uux + hx = -^{u''^)x, ^^g^ 

hf + uhx + hux = —{h'u')x- 

The right-hand sides of the above system are the subgrid-scale quantities that must be param- 
eterized, i.e. it is necessary to find a functional relation that allows one to express these terms 
using only the grid-scale quantities, i.e. 

-l(^)x = fit,x,u^-^\U''^^), -(h^)x = git,x,u('^\U''^^), 

where ri,r2,qi,q2 G Nq. It is the purpose of this section to give two examples of forms of the 
parameterization functions / and g that lead to closed equations 

ut + uux + hx = f, ht + uhx + hux = g (19) 
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possessing certain conservation laws. Here and in what follows, we omit the bars over averaged 
variables in the closed class of shallow-water equations since there is no risk of confusion as only 
averaged forms of dependent variables and their derivatives arise in such equations. 

This notation is also convenient since in several of the examples given below, we consider the 
related problem of finding right-hand sides / and g in system (19) preserving certain CLs of 
system (17), i.e. we consider structure-preserving extensions of the shallow-water equations as 
discussed in Remark 1. 

Before determining ansatzes for / and g that possess CLs or preserve certain CLs holding 
for the free shallow-water equations (18) it is instructive to determine the CLs of system (17). 
We restrict ourselves to multipliers that depend on t, X, U and H, i.e. = A^(t,x,U,H) 
and = A'^{t,x,U, H), where U and H are arbitrary functions of t and x. The determining 
equations (4) for multipliers A^ and A^ in this case become 

Eu{A'Ai{t, X, [/«,//(!)) + A'A2{t, X, i^W)) ^ 0, 
E^^(AiAi(t,x, [/«,/?(!)) + A2A2(t,x,C/W,/fW)) = 0, 

where the Euler operators Ejj and Eh are given by 

Eu =du - Btdut - 'Dxdu^ H , Eh = dn - - B^Oh^ H . 

Splitting the system (20) with respect to the derivatives of U and H, one obtains the following 
system of determining equations for CL multipliers 

A]j-Al = 0, AIj-HA% = 0, a? + [/a2 + A^ = 0, Al + UAi + HAl = 0. (21) 

Differentiating the third equation in system (21) with respect to H and multiplying the resulting 
equation with H, one obtains, upon recombining with the first two equations, the equation 

A^;, + UAl, + HAl, = 0. (22) 

Differentiating the last equation in system (21) with respect to U and then combining the 
resulting equation with (22), one finds that A], = 0. Then differentiating the second equation 
in (21) with respect to x one gets A'j^^ = and thus A'jj^ = in view of the third equation in 
system (21). Differentiation of the fourth equation in (21) with respect to x yields A^^. = and 
hence A^^, = due to the third equation. 

Using these results, the integration of system (21) leads to 

A^ = -citH + X^(U,H), 

(23) 

A^ = ci{x-tU) + X'^{U,H), 

where ci = const and = X^(U,H) and = \^{U,H) are any functions satisfying the 
system \\j — A^ = and \\j — HX^ = 0. Hence there are an infinite number of associated 
CLs, which reflects the possibility of linearizing the quasilinear system (17) using a hodograph 
transformation (interchanging the dependent and independent variables), see also Section 3.3. 
More details on the connection between CLs and linearization of partial differential equations 
can be found in [3, 10]. 

3.1 Conservative parameterization schemes via direct classification 

In this subsection we give an example for the construction of conservative parameterization 
schemes using the technique of direct CL classification. 

As an example, consider the problem of finding diffusion terms of the form 

ut + uux + - F{h,Ux,hx)uxx = 0, ^^^^ 
ht + uhx + hux = 0, 
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allowing for CLs arising from multipliers of the form = A^(i, x, U, H) and A^ = A^(t, x, U, H), 
where U and H are arbitrary functions of t and x. The problem is to first determine correspond- 
ing CLs arising for arbitrary F. Following this, we determine particular forms of F that yield 
additional CLs. The classification is done up to equivalence. 

Theorem 3. The equivalence algebra 0~ of the class of one- dimensional dissipative shallow- 
water equations is generated by the following basis elements 

dt, dx, tdx + du, tdt - udu - 2hdh - Fdp, xd^ + udu + 2hdh + 2FdF- (25) 

Besides the associated continuous equivalence transformations (25), the class of equations (24) 
admits two independent discrete equivalence transformations, which are given by {t, x, u, h, F) i— )• 
{—t,—x,u,h,—F) and {t,x,u,h, F) i— >• {—t,x,—u,h,—F), respectively. The continuous and dis- 
crete equivalence transformations form the equivalence group G~ of the class (24) . 

The determining equations for CL multipliers A^ and A^ are given by 

Ajj-Al = 0, Ah-HAl = 0, Al + UAl + A'^ = 0, 
A^Fh^ = 0, H^A^Fhu^ + iU.A}; + H.Ajj + Al)Fu^ + 2A^F = 0, 
2H,A'Fhh. + 2{UxA\j + HxA\j + AI)Fh^ + A^F^^ + FA]j = 0, 
HxA^Fhh^ + {AI + UxA\j + HxA]j)Fh^-A^Fh-A\jF = Q, 
HIA^Fhh + 2{UxH.,AIj + H^A], + H.,AI)Fh + {U^Al^ + 2UAAuh + 
HlA]jH + 2UxA]j, + 2if,.Ak. + AL)i^ + UAl + HAl + A^ = 0. 

Assuming F to be arbitrary, one can split the above system with respect to the various 
derivatives of -F, which then leads to the solution A^ = and A^ = const. That is, for arbitrary 
F , the only CL admitted by system (24), corresponding to the specified class of multipliers, is 
conservation of mass. 

We now consider particular forms of F for which system (24) possesses additional CLs. From 
the classifying condition A^Fh^ = in system (26) it follows that A^ = when Fh^ / 0. In this 
case, no CL extension exists. Thus, we only study the case of Fh^ = subsequently. In this 
case, the system of determining equations (26) simplifies significantly since all terms involving 
derivatives of Fh^ vanish. Consequently, it is possible to split the resulting equations with 
respect to the powers of H^. The resulting system of determining equations is given by 

A],-Al = 0, AIj-HA% = 0, A^ + [/A2 + a? = 0, (27a) 

A^Fh + FA]j = 0, (27b) 

iU^Ah + Al)Fu^ + 2AIjF = 0, (27c) 

{U^Ahu + 2U,AIj, + Al,)F + UAl + HAl + = 0. (27d) 

Equations (27a) do not involve the constitutive function F and thus can be integrated immedi- 
ately. This results in 

A^ = -^C2tf/^ + (ciH + C2X + C3)U - C2tHlnH - ai{t)H - ^a'l{t)x^ 
-a'2{t)x + a3{t), (28) 
A^ = -ciU^ - {ai{t) + C2t) U - {c2{Ut -x)-C3)lnH + ciH + a[{t)x + a2{t), 

where ai{t), 02 (i) and 03 (t) are arbitrary smooth functions of t, a prime denotes the derivative 
with respect to t and Cj, i = 1, . . . , 3, are arbitrary constants. 
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Each of equations (27b)-(27d) explicitly involves the constitutive function F and these four 
equations are solved using compatibility analysis. Thus, different cases arise. 

Case (I), Fh = 0. From Eq. (27b) it follows that A]^ = and Eq. (27a) imphes that Afj = 0. 
From the form of the multipliers (28) we find that ci = C2 = 0, ai{t) = and thus the multipliers 
in this case are of the form 

= C3U - a2{t)x + asit), = C'slnH + a2{t), 

Substituting this into Eq. (27d) results in 

a2it)U + a2{t)x - a'^it) = 0, 
which leads to 02 (i) = const = 02 and 03 (t) = const = 03. Eq. (27c) then reduces to 

C3iU,Fu^+2F) = 0. 

This equation implies that either (i) C3 = and F = F{Ux) is arbitrary or (ii) C3 7^ and 
F = cq/U"^, with being the integration constant. Since cq 7^ by assumption, we can use the 
transformations from the equivalence group G~ to scale cq = 1. Thus in subcase (i) there is one 
additional CL which is of the form 

^tu + Dx(]-u^ + h- [ F{ux)dux] = 0. (29) 



.2 

In subcase (ii), the multipliers are 

A^ = C3C/ + a3, A^ = C3lnF + a2 

so that here there are two additional CLs. The CL associated with 03 = const is CL (29) 
provided that F = l/u^. The second CL associated with C3 is of the form 

Bt Qii^ + (/iln/i - /i) - + /iln/i+ -^^ =0. 

Case (II), fn / 0. Substituting A^ in the form given in Eq. (28) into Eq. (27b) we can split the 
resulting equation 

^C2tC/^ + {ciH + C2X + cz)U - C2tH\nH - ai{t)H - ]^a'{{t)x'^ - a'2{t)x + az{t)^ Fh 

+ [ciU - C2t(l + \nH)- ai{t)) F = 0, 

with respect to powers of U and x since F = F{H,Ux) only. Splitting with respect to C/^, 
x'^ and X implies that C2 = 0, ai{t) = a\t + and a2{t) = const = 02, respectively, where 
a\,ai = const. Differentiating the simplified equation twice with respect to t leads to a'^{t) = 
or 03(4) = a^jt + 03, a3,a3 = const. The above equation thus simplifies to 

{{ciH + C3)U - {a\t + a\)H + a\t + a^) Fh + {ciU - a\t - a\) F = 0. (30) 

Plugging the simplified form of A^ and A^ into Eq. (27d) leads to a\ = 0. Splitting Eq. (30) 
with respect to t and U yields the following system of three equations 

a\HFH + a\F = 0, {ciH + C3)Fh + ciF = 0, (a?i7 - a^Fn + a^F = 0. (31a) 

At the same time, the remaining classifying equation (27c) gives 

{ciH + C3)(2F + UxFuJ = 0. (31b) 
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System (31) allows one to find four inequivalent solutions. One either has (i) {ci,C3,a\) = 
(0,0,0), (ii) (ci,C3) / (0,0), a} = and a? = 6ci, = -6C3, (iii) (ci,C3,Q!^) = (0,0,0) and 
a5 = —Sa\ or (iv) (03,03) = (0,0), ci 7^ and a\ = 6ci, = eci, 5,e = const. In all other 
cases F = 0, which is excluded from consideration. In subcase (i), the classifiyng equations (31) 
are integrated to give 

9{Ux) 



where g = g{Ux) is an arbitrary non- vanishing smooth function of and 7^ since otherwise 
the assumption Fh 7^ would be contradicted. In this subcase, the multipliers and are 
given by 

The two CLs that correspond to these characteristics are conservation of mass, which follows 
from the multipliers (A^,A^) = (0,02) and 

Dt{aihu — a^u) + D^; ^-a^/i^ — 2*^3^^ ~'~ '^1^^^ ~ cks^ ~ J 5'('"x)dna;^ = 0, 

for the multipliers (A\ A^) = (alH - a'^,a1U). 

In subcase (ii), the solution of (31) is F = cqU^"^ / {ciH + C3) and using the equivalence 
transformations from C one can put cq = 1. The multipliers A^ and A^ are of the form 

A^ = {ciH + C3){U -5), A^ = ^^ciU"^ - 5ciU + c^hiH + ciH + a2. 

The CLs associated with the found multipliers are again conservation of mass associated with 
(Ai,A2) = (0,02) and 

V)t{cihu + C3u) + (^{C3U^ + cih^) + {cm^ + Oi)h + — ) = 
V2 UxJ 

for the multiplier (A^,A^) = {ciH + C3,ciU). The third CL, corresponding to the multipliers 
(Ai,A2) = {{ciH + C3)U,ciU'^ /2 + C3\nH + ciH), is 

Q(ci/i + C3)u2 + C3{h\uh-h) + ]^cih? - + Qci/in^ + ^czu'^ + 

c^huhih + cih?u -\ ) =0. 

UxJ 

In subcase (iii) the solution of system (31) is F = H~^g{Ux), for an arbitrary non- vanishing 
smooth function g{Ux)- Here the multipliers are of the form 

A^ = a\{5- t)H, A^ = a\{x - Ut) + 5a\U + 03. 

The three CLs associated with these multipliers are again conservation of mass (A^ = 0, A^ = 
02), the CL 

Dt{xh — thu) + Da; (^hu — thv? — ^^h^ ^ ^ J ^(^^'^"^^^^ ~ ^' 
corresponding to (A^, A^) = {—tH,x — Ut) and the CL 
Bt{hu) + + hu^ - j g{ux)du^ = 
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which stems from the multiphers (A^, A^) = {H, U). 

In the final subcase (iv) we obtain F = coH~^U~^ from the integration of the system (31). 
Again cq = 1 mod C. The multipliers in this case are 

A^ = {ci{U - 6t) - eci)H, A^ = ^c^U^ + ciH + 6ci{x - Ut) - eciU + aa- 
Besides the obvious conservation of mass, the three other CLs include 

Di(x/i - thu) + (xhu - thu^ - ]-th^ - — ) = 0, 
stemming from (A^, A^) = {—tH, x — Ut), 

Dt{hu) + D^. (Ih^ + hu'^ + —) = 
which is obtained from the multipliers (A^, A^) = {H, U), and conservation of energy 

associated with the multipliers (A^ A^) = {UH, [7^/2 + H) 

Comparing the results of the two cases Fh = and Fh 7^ we have proved the following. 

Theorem 4. For F ^ any equation from the class of dissipative systems of one- dimensional 
shallow-water equations (24) has at most four linearly independent conservation laws arising 
from multipliers of the form = A^{t, x, U, H) and A^ = A^(t, x, U, H). A complete list of G^- 
inequivalent equations and their associated conservation laws is given in Table 1, where g = g{ux) 
is an arbitrary non-vanishing smooth function of Ux and c\ and C2 are arbitrary constants with 

3.2 Conservative parameterization schemes via inverse classification 

In this subsection, two examples are presented that illustrate the procedure for finding conser- 
vative parameterization schemes through inverse group classification. 

Parameterizations conserving energy, mass and momentum. In this example we focus 
on four physical CLs, namely conservation of mass-specific momentum, mass, momentum and 
energy which correspond to the multipliers 



A^ = Cl + c^H + aUH, A^ = C2 + C3C/ + C4 Q^72 + , 



(32) 



for ci,...,C4 G M for the shallow-water equations (17). It can be checked that the above 
multipliers (32) satisfy the system of multiplier determining equations (21) and thus yield CLs 
for the shallow-water equations (17). The canonical forms D^p-l-D^X = of the associated CLs 
are 



Cl : 


P = 


u, 


X = 


+ h 


mass-specific momentum 


C2- 


P = 


h, 


X = 


uh 


mass 


C3: 


P = 


uh, 


X = 


U^h + 


momentum 


C4: 


P = 


^{u^h + h^) 


X = 


^hu^ -\- h'^u 


energy. 
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Table 1. Conservation law classification of a class of one-dimensional dissipative shallow- water equations 
Form of F Conservation laws 

VF CLi = Dt/i-FD^(/iii) = 



p _ g{ux) CLi - 



jcih + C2) Btjcihu + C2u) + Qci/i^ + ^C2U^ + cihu^ + C2-J g^du^^ = 



(cih + C2)ul 



CLi = 0, 

Bt{cihu + C2u) + \ \{c2u'^ + cih?) + (ciu^ + C2)h — ) = 0, 

(^{cih + C2W + C2{h\nh-h)+'^cih'^ -t^ + 

Dj. [ —cihv? -\ — C2V? + C2hulnh + cih?u + — 
V2 3 Ux_ 

CLi = 0, 

Dt{xh — thu) + Dx {xhu — thv? — \th^ + t J g{ux)dux) = 0, 
Bt{hu) + + hu^ - J g{u^)dux) = 

CLi = 0, 

^ 1 Dtixh — thu) +Dx ( xhu — thu"^ — -th"^ — 



hu^ 



T)t{hu) + ( + /in2 + ^ ) = 0, 
V2 UxJ 

Dt [\hu'^ + \h'^ - t ) + [\hu^ + h'^u + — 
V 2 2 / V 2 Ur. 



We now consider the problem of finding parameterization schemes (19) that preserve the 
above four multiphers. For the sake of demonstration we hmit ourselves to constitutive functions 
of the form 

/ = f{x,u,h,Ux,hx), g = g{x,u,h,Ux,hx), 

where here and in what follows we omit the bars and the averaging of the dependent variables 
is to be understood. That is, we look for functions / and g that satisfy 

Ai(Ai - /) + A2(A2 -g) = Btp + D,,Xo, (33) 

for the four multipliers (32) where Ai = ut + uux + hx, A2 = ht + uhx + hux as before. Note 
that adding input terms to the shallow- water equations (17) will in general lead to modified 
expressions for X, which is why we use Xq in the CL (33). However, the conserved quantity 
p remains unchanged since the parameterization functions do not depend explicitly on t or 
derivatives of the unknown functions with respect to t. 

Applying separately the Euler operators and E/^ with respect to U and H of Eq. (33), 
one obtains the system of determining equations for CL multipliers. Since the multipliers A^ 
and A^ are already prescribed, the resulting system 

E[;(Ai(Ai - /) + A2(A2 - g)) = 0, ^h{^\^i - f) + AH^2 - 9)) = 0, 

is now the system of the determining equations for the parameterization functions / and g. The 
determining equations can be split with respect to the constants ci, 02,03,04 and the uncon- 
strained variables which are t, Ut, Ht, Utx, Htx, Uxx and Hxx- Splitting with respect to the 
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highest derivatives arising yields the elementary equations 

which can be integrated to give the following constrained form for the functions / and 

/ = f\x, U, H)a, + fix, U, H)H, + fix, U, H), 

g = g\x, U, H)U^ + g\x, U, H)H., + g'ix, U, H). 
The remaining determining equations can be split to yield the system 

f = g' = Q, f = f = 0, g' = f, 

g'-Hf = 0, f^-f = 0, fl^-HfH-f = 0. 

System (34) can be integrated to give the most general form of the functions / and g admitting 
CLs of interest. In particular, one obtains 

f = fU, + fH,, 
g = fHU, + fH,. 
A particular class of solutions is given by 

f = (oi sin VbU + 02 cos VbU) (asJoilVbH) + 04^0(2%/^)) + 05, 

1 / X (36) 

/ = —J=icil COS VbU — 02 sin VbU) (asJiilVbH) + aiYii2VbH)) . 
\ H ^ ' 



In this solution, oi, . . . , 05, 6 = const, 6 > and J„ and Yn are the Bessel functions of the first 
and second kind, respectively. 

The form (36) for f and f does not lead to a particularly physical parameterization ansatz. 
Physically more relevant forms for / and g can be found upon imposing other restrictions on 
these functions, which involves finding another interesting set of solutions of (34). An example 
of such a construction is the subclass of parameterization schemes of the form (35) which in 
addition to (34) satisfies the equation = 0. In this case, the functions f and f in (35) are 
given by 

f = P,lnH + h. f = ^^^, 
where /3i,/32,/53 are arbitrary constants. 

Conservative momentum dissipation schemes. As discussed in Remark 1, apart from the 
problem of finding conservative parameterization schemes another question of physical impor- 
tance is to construct input terms that preserve some of the geometric structure of the initial 
system. 

We illustrate this idea by constructing the most general dissipation term of the form 

ut + uux + hx = fix, h, hx,Ux,Uxx), ht + uh^ + hu^ = 0, 

for the shallow- water equations that preserves the multipliers A^ = cq and A^ = ci, i.e. we set 
/ = fix, h, hx,Ux,Uxx) and g = in system (19). We do not aim at conserving momentum or 
energy in this case, because our aim in this example is to construct a dissipation for the shallow- 
water equations, which by definition should violate energy and momentum conservation. 

Using the same procedure as outlined in the previous example, we find that / should be of 
the form 



/ = f'uxx + (y fldux + f^^ ^- + 1 [I f^hdux + dh + CUx + f 



(37) 



where c G M, = f^ix,h,Ux), = f'^ix,h) and = f'\x). 

From the form of (37), one observes that the requirement of conserving both mass and mass- 
specific momentum leads to a quasi-linear dissipation scheme, i.e. / is linear in terms of Ux 



''XX • 
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3.3 Conservative and invariant parameterization 

We now turn to the problem of finding parameterization schemes that are both conservative and 
preserve certain symmetries of the original (unaveraged) system of differential equations. We 
illustrate this idea with the two examples of the previous section. 

The maximal Lie invariance algebra g of the system of one-dimensional shallow-water equa- 
tions (17) is infinite dimensional and has the following basis elements 

tdt + xdx-, xdx + udu + 2hdh-, td^ + du, 

n n (38) 

{2x - Qtu)dt + {Qh - 2,u^)td:, + {u^ + Ah)du + '^hudh, T{h,u)dt + C{h,u)da^, 
where the functions r and ^ run through the set of solutions of the system 

Ch - UTh + Tu = 0, Cu- UTu + hTh = 0. 

The infinite dimensional part of g indicates the existence of a linearization transformation, which 
is given through the hodograph transformation, i.e. t = t and = x are the new dependent 
variables and u and h are the new independent variables. 

The question of which symmetries one aims to preserve when constructing a (conservative) 
parameterization scheme should be answered using physical arguments. For example, processes 
that are to be parameterized in the framework of classical mechanics should be represented in 
such a manner so as to be invariant with respect to the Galilean group. Choices for subgroups 
to be preserved by a parameterization scheme can be also motivated from compatibility with 
certain boundary- value problems, see also the related discussion in [7]. 

As outlined in Section 3.3, when constructing parameterization schemes that are required 
to be both invariant and conservative one can follow two ways, namely using direct or inverse 
group classification. We now use the examples worked out in the previous section to illustrate 
both ways. 

Invariant conservative parameterization using direct classification. The direct group 
classification method can be illustrated with the first example from the previous section. Essen- 
tially, the only freedom one has left with the parameterization (35)~(36) is to set to zero some 
of the constants ai, . . . , as and investigate which symmetries the resulting systems have. This 
analysis should be done up to equivalence of the class of equations of the form (35)-(36). 

The kernel of the maximal Lie invariance algebra from this class is given by the subalgebra 
of consisting of the basis elements 

0^ = {tdt + xdx,dt,da:). 

It turns out that the only extension of this kernel algebra arises when ai = 02 = as = 04 = 
and 05 = c 7^ 0. However, in this case the transformation u = u — c maps the resulting case to the 
initial shallow-water equations (17), i.e. the parameterization becomes trivial. Stated in another 
way, the constant 0:5 is inessential for the classification problem of the class (35)-(36), which 
has the two inequivalent solutions (01,02,03,04) 7^ (0,0,0,0) and (01,02,03,04) = (0,0,0,0). 
Only the first solution leads to a nontrivial parameterization scheme. 

A similar analysis could be carried out with the other parameterization schemes of the 
form (35) that can be constructed by finding only particular solutions of the system (34). As 
in the previous case, these classes of parameterizations essentially only depend on certain con- 
stants. The classification problem then basically reduces to finding those constants that can 
be set to zero by a proper transformation of the equation variables (inessential constants) and 
studying the classification problem with respect to the remaining (essential) constants. This is 
a straightforward task and is not considered further in this paper. 

Invariant conservative parameterization using inverse classification. We now focus 
on the problem of finding elements of the class of shallow- water equations (19) with / given 
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by (37) and g = that are invariant under a certain subgroup of the maximal Lie invariance 
(pseudo)group G of the shallow-water equations. That is, we use the inverse symmetry classifi- 
cation strategy to find models from this class. As a symmetry subgroup G^, we single out the 
four-parameter subgroup of G that is generated by the four-dimensional Lie subalgebra of q 
with the basis elements 

dt, dx, tdx + du, tdt + {I + d)xdx + dudu + 2dhdh, 

for an arbitrary constant d G M. The reason for choosing this particular subalgebra is that 
for physical arguments we require our dissipation scheme to be invariant under the Galilean 
group (generated by the first three elements of g^) and to have a scaling symmetry. We only 
require invariance under a single scaling instead of the two scalings admitted by the original 
shallow- water equations (17) as adding dissipation to a hydrodynamical system usually breaks 
one scaling symmetry. For example, the shallow-water system with classical linear dissipation, 
/ = Uxx has the above scaling symmetry provided d = —1/2. 

As discussed in Section 3.3 one can use the replacement theorem to construct parameterization 
schemes that are both conservative and invariant. To this end, we determine the second order 
differential invariants of the Lie algebra g^. They can be found by using infinitesimal techniques 
through prolongation of to the action on second derivatives of u and h, and then invoking 
the infinitesimal invariance criterion [8,10,11,21,23]. These invariants can also be found using 
the moving frame method [12, 15, 16]. In particular, on the space of equation variables for the 
class (19) with / satisfying (37) and (7 = 0, which is the subspace of the second jet space 
with coordinates {t,x,u,h,ut,Ux,ht,hx,Uxx), there are five elementary invariants given by 

h = hul'', h = hxut\ h = UxxUx^^+''\ ^^^^ 

h = ui~^{ut + uux), h = uif'^ih + uhx). 
The invariant representation of system (17) is therefore 

U + h = Ut^^l, h + h= 

from which it follows that the multipliers FJ, F^, Y:\ and F| in system (9) are F} = n^-^, V\ = 0, 
F^ = and Tl = uf-^. 

The system of shallow-water equations with the dissipation scheme (37) is invariant under 
the group generated by the elements of provided that 

Fi/ = /(/i,/2,/3) (40) 

holds for some function / that can depend at most on the invariants /i, I2 and Is since none of 
the functions and in (37) depends on ut or ht. 

We now determine some dissipation schemes that fulfill the above requirement. For the sake 
of simplicity, we assume that = = 0, since we are mainly interested in the second-order 
term proportional to Uxx and thus in finding functions that lead to invariant and conservative 
diffusion schemes. As none of the invariants /i, I2 and I3 depends on x we have that = 0. 
The function / should depend linearly on to match with Uxx in /• Comparing the coefficients 
of Uxx we find that 

Note that a cannot depend on I2, since does not depend on hx- The remaining condition 
that has to hold is that 

hx [ fldUx=ul-''p{h,h). 



22 



The function (3 now cannot depend on as there is no u^x term in the above left-hand side. 
Since the left-hand side is linear in /i^., we find that /3 = /27(A) and thus get 



which imposes a relation between the functions a and 7. This relation is particularly straight- 
forward to evaluate for polynomial functions a. To give an example, let us set a = 2cdlf, for 
c G M. This leads to the shallow-water equations with dissipation in the form 



The usual linear dissipation / = I'Uxx, G M, falls into this class when putting d = —1/2 
and using a(/i) = u. 

4 Conclusion 

In this paper we have studied the problem of finding physical parameterization schemes that lead 
to closed systems of averaged differential equations which possess nontrivial local conservation 
laws. A main motivation for our work is that by its formulation, one cannot expected to find an 
exact solution to the parameterization problem. The determination of the entire subgrid-scale 
structure of a real-world process when one has at their disposal only the grid-scale information 
is not feasible for nontrivial physical problems. Hence, any auxiliary information that can be 
used to limit the possible form of a parameterization scheme by imposing some physically and 
geometrically relevant structural constraints is highly useful. The preservation of symmetries 
and conservation laws can serve as such relevant constraints since these two properties are closely 
linked with the physics encoded in a system of differential equations. 

A systematic toolbox of methods that allows one to systematically find parameterization 
schemes with symmetry properties using group classification techniques was formulated in [6,24]. 
As far as we know, this paper is the first to use the analog toolbox for finding parameterization 
schemes preserving conservation laws. 

The results of this paper illustrate that the requirement of preserving a particular set of con- 
servation laws when constructing physical parameterization schemes can lead to rather specific 
forms for these schemes. This is in striking contrast to the case of invariant parameteriza- 
tion schemes, since here there are in general an infinite number of possibilities to construct a 
subgrid-scale closure possessing a prescribed maximal Lie invariance group. Moreover, as to be 
expected, the more conservation laws one aims to conserve when constructing a subgrid-scale 
closure or any other additional model for a system of differential equations, the less freedom 
one has to adjust the closure by including other desirable properties. Thus, the requirement 
of preserving certain conservation laws (and symmetries) can lead to rather specific parameter- 
ization ansatzes which in consequence could potentially simplify the construction and testing 
procedures for subgrid-scale closure models. 

In the present paper, the primary focus of the presentation was to give a careful exposition of 
the different ideologies for finding conservative closure schemes. The system of one-dimensional 
shallow-water equations served as a proof-of-the-concept example but did not reveal new physical 
insights. More realistic examples of conservative parameterization schemes for the governing 
equations of hydro-thermodynamics will be presented elsewhere. 
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